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vM Abstract. The control of condensed matter systems out of equilibrium by laser 

pulses allowrs us to investigate the system trajectories through symmetry-breaking 
' I ' phase transitions. Thus the evolution of both collective modes and single particle 

excitations can be followed through diverse phase transitions with femtosecond 
resolution. Here we present experimental observations of the order parameter 
trajectory in the normal— ^-superconductor transition and charge-density wave 
^i^ ordering transitions. Of particular interest is the coherent evolution of topological 

fli defects forming during the transition via the Kibble-Zurek mechanism, which 

e appears to be measurable in optical pump probe experiments. Experiments on 

CDW systems reveal some new phenomena, such as coherent oscillations of the 
-4_> order parameter, the creation and emission of dispersive amplitudon modes upon 

C^ the annihilation of topological defects, and mixing with weakly coupled finite- 

C frequency (massive) bosons. 
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Coherent topological defect dynamics and collective modes in superconductors and electronic crystals.2 

1. Introduction. 

The idea that non equiUbriuni phenomena in natural systems can be described in terms 
of a temporal evolution of a non-linear Landau expansion of the free energy has lead 
to theories as wide afield as elementary particles to cosmology. In condensed matter 
systems, where the idea originates, the opportunities for the study of the behaviour 
of nonlinear systems abound. However, typically, the study of these systems has been 
limited to near-equilibrium situations, where the system evolves slowly through the 
transition. Yet both particle physics experiments and cosmology distinctly deal with 
the temporal evolution of highly non-equilibrium systems. Typically one detects the 
decay products well after the decay itself, in the aftermath of the symmetry-breaking 
transition (SBT). Collisions of elementary particles, the subsequent creation of a high- 
symmetry intermediate state and its decay via symmetry-breaking interactions takes 
place on timescales which are well beyond the resolution of our current technology. 
Cosmological experiments being out of reach, experimental cosmological studies of 
the aftermath of the Big Bang reveal very complex behaviour in its aftermath, still 
beyond current theoretical understanding. As a result we have little insight into the 
non-equilibrium conditions close to the critical time tc of the SBT when the creation 
of the new state is taking place. 

In condensed matter systems the intrinsic dynamics of elementary excitations 
occur on timescales which are becoming accessible with current femtosecond laser 
technology. Particularly the normal-to-superconducting state transition and charge- 
density wave transition are of fundamental interest, as examples with different 
symmetry properties. Thus, typical single particle (quasiparticle) relaxation times 
in high-temperature superconductors are on timescales of 10~^^ s [T], while electronic 
energy relaxation occurs on timescales of 10""'^^ — 10~"'^^s [5]. Many body collective 
states have similar characteristic timescales. The so-called Ginzburg Landau (GL) 
time TGL, [3] which enters into the time-dependent GL equation is tqi^ ~ l/A^ 
where As is the superconducting gap. For high-temperature superconductors, typical 
estimates give tgl ~ 10""'^'^ — 10~^^ s. In superconducting systems, the collective 
mode is overdamped, which means that it is difficult to discern the relaxation of the 
collective mode from single particle excitations. On the other hand, the characteristic 
timescale of collective mode in charge density wave systems is typically around 
2tt/ujam — 0.5 X 10~^^ s (H [S]. Laser pulses can currently be created with sub 
femtosecond duration, so ample resolution is available for the study of the coherent 
evolution of collective modes in superconducting and charge-density-wave (CDW) 
systems. Moreover, by studying the concurrent evolution of the single particle 
spectrum, it may be possible to investigate the transient state of the underlying 
vacuum. 

Following the suggestion by Zurek|6| of laboratory experiments to test Kibble's 
cosmological model [7], "system quench" experiments were performed in a number of 
different systems, including superconductors |H1[S], Bose-Einstein (BE) condensates 
within a trapped atomic gas [TU] and polariton BE condensates [TT] . Experiments on 
the Kibble-Zurek (KZ) mechanism so far have concentrated on the statistical analysis 
of topological defects left behind by the SBT, on timescales long compared to the 
intrinsic GL time tql. Quench rates in these systems were relatively slow, the fastest 
being around 10*K/s [9 . 

In standard optical Pump-probe (P-p) techniques (including THz probe) the 
response is related to the modulation of the dielectric constant (reffectivity or 
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Figure 1. A schematic diagram of the three-pulse sequence and evolution of 
the superconducting order in a superconductor. The different time delays are 
indicated. The ratio of intensities of the D (blue ball), P(red ball) and p{green 
ball) laser pulses are typically 500:10:1 respectively. The electronic temperature 
decays into the sample as a result of the inhomogeneous excitation by laser light. 
The position of a normal/superconducting (N/S) state front, given by Tg = T^ 
is rapidly moving towards the surface, leaving vortices in its wake. In a CDW 
system, the D pulse forces the system into the high-symmetry state, whereupon 
it undergoes a SBT as it cools. In the process, topological defects are created as 
a result of the inhomogeneity of the laser field. 



absorbance). which can have a number of contributions in electronic systems such 
as superconductors and CDWs which are in one way or another related to the 
order parameter ■0: (1) the response due to hot carrier energy relaxation, (2) 
the quasiparticle (QP) recombination across the SC or CDW gap and (3) QP 
recombination across a pseudo-gap which is often present in these systems and 
(4) coherent phonon oscillations due to displacive excitation or impulsive Raman 
excitation. The coherent phonon oscillations cannot be distinguished from the 
oscillations of the order parameter with P-p spectroscopy. Moreover, different phonon 
modes are coupled to ip to various degrees, and the dephasing is different for each 
mode, thus adding further to the complexity of the response from which it is difficult to 
extract the order parameter. We show that using a 3-pulse technique, these responses 
can be isolated to some extent, allowing us to study the coherent evolution of the order 
parameter with time though an SBT with very high temporal resolution. The method 
allows us to study the interactions of the collective mode of the CDW with other 
modes of different symmetry as well as annihilation of topological defects, revealing 
new finite frequency dispersive Higgs-like field excitations released as a result of domain 
wall annihilation. 

2. Quench experiments with a multiple pulse all-optical technique. 

The principle behind the technique involves the use of a multiple pulse trains to control 
the system and measure the time-response of the order parameter. A schematic 
diagram of the pulse sequence is shown in Fig. 1. A strong laser destruction (D) 
pulse - whose intensity is carefully adjusted to cause a perturbation of appropriate 
magnitude - is used to rapidly transfer the system from an ordered (broken symmetry) 
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Figure 2. The calculated evolution of the order parameter V-d{*) through the 
transition after an intense destruction D pulse at i = (dashed blue line). When 
an additional P pulse is present, the trajectory tpopii) (solid blue line) can be 
drastically modified, depending critically on Ati_2. The difference between the 
two is highlighted by the shaded areas. The evolution is calculated for three 
different time delays Aii_2 = 350, 400 and 450 fs. The intensity of the P pulse is 
10% of the intensity of the D pulse. The calculation in this figure is done with an 
exponential perturbation function fi{t) for the D pulse (shown by the dashed red 
line) and fi0p{t) with both pulses (solid red line). The green line is the response 
function in reflectivity predicted according to Eq. |10| 



state to the disordered (high-symmetry) state, whereafter the system reverts back to 
equilibrium through the SBT. The state of the system at any given time after the D 
pulse is determined by standard pump-probe spectroscopy sequence, which involves 
first exciting the system by an additional weak perturbing P pulse delayed by Aii_2 
with respect to the D pulse, and subsequently measuring the resulting change of 
reflectivity AR/R of the sample by a weak probe (p) pulse (see Fig. [ij. In sec. 3.1 
we will show how the optical response AR/R can be related to the order parameter 
in superconductors and CDW systems. 

On the level of field theory, the trajectory of the system through the SBT 
can be modeled within time-dependent Ginzburg Landau theory with a free energy 
functional corresponding to the symmetry of the problem in hand. For a charge 
density wave system, to capture the salient physics, we can assume an order parameter 
\I/ = A(i,r)e"^ and make the assumption that the relaxation on phase (j) is slow 
compared to tqi^, which means that the trajectory is dominated by the dynamics of 
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A(t,r) da [13]: 

d^A dA d^A 

where ujam is the angular frequency of the collective amplitude mode, ^ is 
the coherence length, F = Aojam/^am is the damping of the AM, and a(t, r) is 
temporally and spatially dependent and is derived from the usual GL temperature 
a(i,r) = 1 — l^{t,^) describes the control parameter. fJ,{t,r) may be T(t,r)/Tc, as 
given by original GL theory, or may be approximated by an exponential function 
which signifies the cooling of the system. The spatial variation of the light intensity 
is accounted for by an excitation function ^{t, z) = 'j^'^' exp(— z/A) where A is the 
optical penetration depth. Using experimental values for vam = '-^am/'^tt = 2.4 THz, 
the line-width Ai^am = 0-05 THz and coherence length ^ = 1 nm [T^ we can compute 
A{t,z). 

In 3-pulse experiments, the P pulse presents an additional perturbation which 
can modify the trajectory, especially if it occurs close to the critical time tc- In Fig. 
l2]we show calculated trajectories of ^p for three different delay times Aii_2 of the P 
pulse with respect to the D pulse. The "butterfly effect" is very evident: V msiy end 
up as either +1 or -1, critically depending on Aii_2. 

The equation of motion for a superconductor within the time-dependent 
Ginzburg-Landau (TDGL) model is more complicated because the excitations are 
charged and ip is coupled to the electro-magnetic fleld, and are usually considered to 
be overdamped, so they do not include any second-order time-derivative: 



, i$^ ) = -ar(i,r)V'-V'IV'r- (»V + a)^i/^ + ry (2) 



V^$ 



*(V'*VV'-V'VV'*)+a|i/'|^ . (3) 



.2 

These equations have been used for analytic and numeric analyses of dynamics of 
topological defects (vortices and phase slips) in superconducting wires, as reviewed in 
|14[[T5]. The vector potential a is written in units of ^^ {cpo is the flux quantum, ^ is 

the coherence length) and the electrostatic potential $ in units of 2^, with e being 
the elementary charge, Tg is the transverse relaxation time and h the reduced Planck 
constant. The only dimensionless parameter left in the equation is the ratio u — -^ 
between the longitudinal relaxation time and the transverse relaxation time. Typically, 
for high-temperature superconductors rg « 1 ps and u ~ 5 |15j . The Langevin noise 
term rj introduces microscopic fluctuations into the dynamics |15j . 

3. Control of the quench rate and optical response function 

In the simplest approximation, suggested by the usual a parameter in the GL 
equations, the excitation function /i can be the system temperature. In the case 
where the electronic system is being discussed, such as superconductors and CDWs, 
it is the electronic temperature T^it) which is the relevant control parameter. In 
recent years, significant progress has been made in understanding the time evolution of 
electron and lattice temperatures in superconductors and CDWs following excitation 
by a laser pulse. The quench rate is determined by the energy relaxation rate of 
electrons at early times up to ~1 ps, and the lattice Ti(i) on longer timescales (tens 
of picoseconds and beyond). The time-evolution of Te(i) and Ti^{t) respectively, which 
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govern the quench process, may be estimated using a two-temperature model |16| : 
leTe^ = -"fL{T,^n)+Eit) and C^T)'^ = TlITo-Tl), where E{t) is a function 
(Gaussian) describing the energy density per unit time suppHed by the laser pulse. The 
pulse energy Eu = Epjd^ where Ep is the pulse energy. We used the temperature- 
dependent thermal constants jeiT) and Cl(T) from experimental datajTS] and the 
measured energy loss rate 7l — 340 K/ps for Lai.gSro.iCu04 [2]. The dependence of 
Tc{t) and Ti^{t) on the energy density per pulse U — J E{t)dt is shown in Fig. ^) for 
the first few picoseconds. The red contour lines show when either T^ or Tl crosses T^. 
Initially, the rate of cooling 7e — ^^ is very rapid until Tg reaches T^. Thereafter, 
Te follows Tl, and is dominated by phonon anharmonic decay or phonon escape and 
thermal diffusion processes, all of which are much slower than electronic energy loss. 
Thus, to change the quench rate we can adjust the energy density Ud of the D pulses, 
so that the system either cools rapidly through electron thermalization, or more slowly 
by phonon decay and diffusion processes. As shown in Fig. Isk), for low Ud, only Tg 
exceeds Tc, so the quench back through Tc will be fast and dominated by electronic 
cooling. For larger Ud, Tl exceeds Tc, so on cooling the quench rate though Tc is 
significantly slower and dominated by lattice cooling. 

3.1. The response function 

The key issue is detection of the order parameter through the SBT. Here we show 
how for superconductors and CDW systems, the response measured in 3-pulse 
experiments is related to the order parameter. While a full kinetic theory which 
would take into account single particle and phonon relaxation processes through 
the transition is beyond reach at present, we can obtain good approximations for 
the quasiparticle response functions in superconductors and collective mode in CDW 
systems respectively based on well-accepted phenomenological theory. 

The temperature dependence of the difference of the optical reflectivity between 
the superconducting state R^ and the normal state i?„ in a superconductor has been 
found to be described very well by an expression derived from the Mattis-Bardeen 
formula J20]: 



A(T) 



2 



RsiT)-Ra^A \' In 



1.47a; 

XiT) 



(4) 



where A(T) = 1 — (T/Tc) is usually used to describe the temperature dependence 
of the gap, bj is the frequency of light in units of h, Rn is the reflectivity in the normal 
state. We need to calculate the transient change of reflectivity 6R, as the system 
is evolving through the transition in time, so we explicitly replace A(T) with A(t). 
Using the 2-fluid model, we substitute A^ — AqUs — Aq(1 — rig), where rig is the 
superfluid density and riq is the quasiparticle density, and Ao is the gap at T = 0. We 
can relate SR to the photo-excited carrier density np using the fact that Sn 
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1.47cj 



A it) 

If we ignore the derivative of the logarithmic correction with respect to A 
substituting Aq for A(t) in the logarithm, we obtain Rs{t) — i?„ ~ A—^ 
6R then simplifies to: 

,„ dRdA^ 2AA{t)^^ AAl 

oA on„ w^ (jj^ 
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Under near-bottleneck conditions, when the QPs and the phonons are in near- 
equilibrium. Up is given by |21| : 

^ i/(A(r(t))) + r(t)/2) 

" 1 + i?v/2T(i)/(A(T(t))Exp[-A(T(f))/r(t)] ' 

where A is the superconducting gap, B ~ v/N{Q)Klc, where v is the effective 
number of phonon modes of frequency Mlc per unit cell participating in the relaxation 
process. Here we have explicitly written the temperature T to be time dependent. B 
can be determined by fitting the temperature dependence of AR/R and N{0) is the 
density of electronic states at the Fermi energy. Substituting Aoip{t) for A(t) from 
the solution of Eq. [3) we obtain: 



sR/Ro.^(i-2iJ'-''-^^ mAo\m\+m/2) 



2^2 



A(t) 



1 + B^2T{t)/Aomt)\Exp[-Aomt)\/T{t)] 



(8) 



The response functions using the appropriate constants for Lai.9Sro.iCu04 are 
plotted in Fig. m together with the time evolution of T'e(i) for E = 1.6 J/cm^, and 
the order parameter \ip{t)\ and |V'(i)P for the case of a homogeneous superconductor 
calculated using Eq. [2] in the absence of fieled a = 0, and $ = 0. The experimental 
reflectivity response SR/R is close to the square of the order parameter \ip\'^ over a 
large range of t, particularly for large ip. The calculated case for an inhomogeneous 
superconductor which takes into account the depth profile of the laser beam is 
discussed in the next section where it is compared with experimental data. 

In CDWs, when discussing the transient response of collective amplitude mode, 
the response is related to the displacement AR of the atoms from equilibrium position 
i?0) where AR ex ijj — 'ijjo, where ^o is the equilibrium value of the order parameter. To 
obtain the optical response, we can expand the dielectric constant near CDW phase 
transition in powers of the order parameter |17| : 

e = eo+C2|^p. (9) 

Here eo is the dielectric constant of the high temperature symmetric phase, C2 is a real 
constant. 

The relevant response in D-P-p experiments can be expressed in terms of the 
difference between the response with and without the P pulse, proportional to Aj^p, 
so to first order |T5]: 

-—^ — ^^lp{t,r,Ati2)~ijlit,r). (10) 

H e 

This response has been tested for the case of the CDW transition in TbTca jT^] . 
We see that by careful design of the experimental probe and the use of 3 pulse 
techniques, it is possible to identify and even isolate the dominant contribution to 
the optical response which is related to the order parameter, thus providing valuable 
information of the temporal evolution of ^ by measurement of either single particle 
and collective excitations through the SBT. 

4. Superconductors: probing vortex dynamics on the picosecond 
timescale. 



In Figure Hk) we show the time-evolution of SR/R in Lai.9Sro.iCu04 measured with 
a. P — p sequence at a sample base temperature of 4 K after a 50 fs laser pulse. All 
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Figure 3. a) The calculated evolution of Te, T^ as a function of D pulse energy, 
b) The order parameter \ip\ and reflectivity response AR/R as a function of 
time through the quench. The calculation of T^ and Tj^ is performed using the 
two-temperature model 1181 1161 using the appropriate specific heat and thermal 
conductivity functions for Lai.9Sro.iCu04 from experimental data [201 . The 
response functions compare the response calculated with Eq. ^ and [6] with \ip\ 
and IV'p. 



laser pulses have a wavelength of 800 nm. The D pulse energy Ed — 0.8 J/cm^ is 
adjusted to heat the sample above Tc- (The threshold values for the destruction of 
the SC state were previously reported in [20J. 

According to Eq. W\ the QP amplitude in P-p experiments are related to |-0|, which 
allows us to compare the data with the calculated trajectory based on a solution of 
Eqs. |2] and [3j The numerical solution calculated without the fluctuation term rj is 
shown by the solid line in Fig. l4k). Comparing with the measured data, we see that the 
prediction is remarkably good overall, except for a depression of the order parameter 
of approximately 10-20 % between 10 and 30 ps. This discrepancy may be attributed 
to vortices which form as a result of the fast quench through T^. ijj ^ in the vortex 
cores, so the spatially averaged value of ip will be depressed if they are present. To 
try and model the effect of vortices we compare calculations with and without the 
fluctuation term 77, while keeping all other parameters unchanged. The introduction 
of fluctuations gives rise to a depression of |-!/)| in the region around 10 ps, shown by the 
dashed line in Fig. HK). Qualitatively, we see that the presence of vortices may explain 
the depression of the order parameter around 10 ps in Lai.9Sro.iCu04. Considering 
the simplicity of both the TDGL model and the crudeness of the approximations 
for the response function (Eq. [?]), the agreement between data and theory is quite 
reasonable. 



5. Topological defects in CDWs: coherent dynamics. 

Transition metal tri-tellurides such as DyTes are excellent systems for femtosecond 
measurements of system trajectories through SBTs to a CDW state. Recently, the 
evolution of single particle and collective excitations though the CDW transition have 
been reported in TbTea [T^]. Here we show the temporal evolution of the order 
parameter in a related material, DyTes, which undergoes a transition to a CDW state 
at 305K. The temporal evolution of the collective mode spectrum after the D pulse 
is shown in Fig. IS] The spectrum was obtained from the Fourier transform of the 
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Figure 4. The recovery of the order parameter ^ after an ultrafast quench 
in Lai.gSro.iCu04. The data points (left panel) represent the dependence of 
the maximum AR/R response as a function of Aii_2, while the solid curve 
is the trajectory of |i/'(t)P calculated from Eq. bj normalized to the value at 
long times after the D pulse. The right panel shows the calculated trajectory 
when fluctuations are included in Eq. [s] via the rj term (dashed line). The 
trajectory without noise r] is included for reference (black line) for the same 
parameters. The difference between the two is attributed to the effect of vortices 
spontaneously created via the KZ mechanism, which annihilate on a timescale 
of 20 ps according to theory predictions. A comparison with the data suggests 
the formation of vortices in Lai.gSro.iCu04 created in the quench on a timescale 
between 10 and 30 ps. Note that we have limited ourselves to the time interval 
where »/) — >■ ipequUibrium^ SO AR/R oc \tp\ is a valid approximation (see Fig. [Sp)). 



transient reflectivity AR/R recorded at each Aii_2 as a function of Ai2_3. We identify 
two modes in the spectra at ^1.68 THz and ^ 2.2 THz at long times Aii_2 > 10 
ps. The latter is assigned to the amplitude mode (AM) of the CDW based on its 
distinct temperature dependence [5^, while the former is a phonon mode (PM). At 
t ~ 1.3 ps. the AM exhibits mixing with the PM. Oscillations of the order parameter 
are clearly visible at short times and appear as modulations of the intensity of both 
modes occurring up to ~ 8 ps. At longer times, the intensity oscillations and the 
frequency recovers to 2.2 THz after 4 x 10® ps. No further relaxation takes place 
beyond this timescale. 

In Fig. lola) we show the calculated 'ip{z, t) from Eq. fllusing parameters for DyTes: 
WAM = 2.2 THz, the AM linewidth F = 0.3 THz, and coherence length ^ = 1.3 nm 
from |24| . For simplicity, the driving term a(t, z) in Eq. nlis assumed to be exponential 
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Figure 5. The temporal evolution of the collective mode spectrum in DyTea 
through the quench into the ordered phase. The plot shows the amplitude of the 
spectrum obtained from the transient reflectivity AR/R as a function of At2-3 
for different Ati_2. The AM shows some intensity around 0.5 THz at short times, 
then crosses a phonon at 1.68 THz at around 2 ps and eventually settles down 
to 2.2 THz at long times. The dotted line shows the recovery of the electronic 
single particle gap determined from the recovery of the SP peak scaled onto the 
u} — Ati—2 plot. At early times, coherent oscillations of the order parameter cause 
fluctuations of the intensity of the AM and the 1.68 THz phonon. Interference 
between the phonon and the AM is clearly observed, including a renormalization 
of the phonon frequency from ~ 1.85 THz at early times (in the high-symmetry 
state) to 1.68 THz in the ordered state. The insert shows the dynamics at short 
times. Distortions of the line-shapes in time are attributed to the presence of 
topological defects created in the quench to the CDW state (see arrow) on the 
basis of comparisons between theory and experiment (see Fig. [ml. The sample 
base temperature is 4 K throughout the experiment. The color scale is normalized, 
so that the maximum amplitude is 1 (red). (In the insert, the scale is expanded 
to show the additional spectral features shown by the arrow.) 



in time and in space: a = exp{—t/T — z/X), where the penetration depth is given by 
the experimental value A = 20 nm [25 . The calculation predicts the formation of 
three domains created within the first picosecond, parallel to the surface, with domain 
walls at ^ 30 nm and ^ 45 nm. After 5 ps, the domains walls are annihilated, leaving 
behind a single domain. The event is accompanied by the emission of a finite frequency 
dispersive amplitude waves traveling towards the surface and into the sample. The 
velocity of the amplitude wave is seen from Fig. [6] to be approximately va — 10 
nm/ps. The wave reaches the surface with 6-8 ps causing a disturbance of the AM 
which is visible as a temporal distortion of the AM spectral line-shape. (The insert 
to Fig. [5] highlights the distortions of cj — Ati2 spectra) . The calculated spectrum 



Coherent topological defect dynamics and collective modes in superconductors and electronic crystals.ll 
100 -, 




Time At^^ [Ps] 



Figure 6. a) The calculated order parameter (colour scale) as a function of 
time and depth into the sample after a laser pulse as a solution to Eq. IT] and 
parameters for DyTes. b) The intensity oscillations are calculated from Eq. |10[ 
taking into account the appropriate boundary conditions and inhomogeneity of 
the laser field in Eq. [l] Distortions to the line-shapes are attributed to the 
annihilation of topological defects, similarly as observed in TbTea 1121 . The color 
scale is normalized, so that the maximum amplitude is 1 (red) and minimum is -1 
(blue). (In the insert, the color scale is expanded to show the additional spectral 
features) 



(using Eq. 10 1 is shown in Fig. b). The predicted distortions have a great deal of 
similarity with the experimentally observed evolution of the spectrum shown in Fig. 
|5] Unfortunately the presence of the phonon interference in this material complicates 
the detailed comparison between theory and experiment. Similar, although more 
pronounced distortions were observed in TbTcs, where the interference from other 
phonons is less problematic |12j . 

In [2(3 it is further shown that incoherent intrinsic topological defect dynamics in 
the related system TbTea occurs on a timescale of ^ 30 ps, similar to the timescale of 
vortex annihilation dynamics inferred from the experiments on Lai.9Sro.iCu04. The 
experiments show that pinned-defect-related annihilation appears to be present on 
much longer timescales of 10^^*^ — 10^^ s. 



6. Discussion 

We have demonstrated that femtosecond optical experiments with multiple pulse 
techniques open the way to detailed studies of the temporal dynamics of systems 
undergoing SBTs and can be interpreted in the context of a phenomenological field 
theory without resort to microscopic theory. The basic analogy with cosmology 
and elementary particle collisions comes from the fact that the SBT take place in 
temporally evolving systems with high temperature intial conditions. Single particle 
fermionic and collective bosonic excitations can be unambiguously identified and 
related to the temporal evolution of the order parameter through the transition and 
topological defects (domain walls and vortices) are shown to lead to experimentally 
observable phenomena. It appears that CDWs are eminently more suitable for 
investigation of topological defects than superconductors because the collective mode 
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conveys significantly more information of tlie system trajectory than single-particle 
excitations, particularly relating to coherent defect dynamics. We have shown that 
by adjusting the laser excitation energy density, the quench rate can be varied and 
controlled, leading to the possibility of investigating domain wall recombination. 
Further studies as a function of quench rate and temperature may be expected to reveal 
systematic behaviour, which can be related to the predictions of the KZ mechanism 
for the generation of topological defects. In the context of the discussion regarding the 
possible cosmological importance of creation and annihilation of cosmic domains and 
strings f7], the present paper demonstrates that coherent creation and annihilation 
of topological defects leads to observable consequences in condensed matter systems. 
Contrary to previous experiments, the distinction between coherent and incoherent 
recombination dynamics, as well as intrinsic and extrinsic defect-related annihilation 
is quite clear, because the rates of decay are vastly different: coherent domain 
recombination occurs on a Tcoh — 3 ps timescale, while intrinsic incoherent domain 
wall recombination timescale is r^ ~ 30 ps. Finally, extrinsic defect-related dynamics 
occurs on timescales of r^xt ~ 100 ps or longer |25]. The experiments thus establish a 
hierarchy of timescales associated with the annihilation of topological defects. 

The evolution of the single particle response as a function of time after destruction 
by a laser pulse in the Lai.gSro.iCu04 superconductor suggests that a depression of the 
order parameter compared with the theory prediction can be understood if vortices are 
included in the within an inhomogeneous model, suggesting intrinsic vortex creation 
and annihilation in Lai.9Sro.iCu04 takes place on a timescale of 10-20 ps. Although 
one might argue that the discrepancy between the trajectory of |'0| and the measured 
response may be attributed to breakdown of the approximations made in deriving 
the response function Eq. [7] these problems are not expected to be important in the 
region around 10-20 ps, but rather at shorter times, closer to tc'- when ip —>■ 1, the 
gap is fully developed, and the QP response is expected to be described well by the 
Rothwarf- Taylor model in the bottleneck regime[23], where the gap A(i) -^ 1 and is 
almost constant, so Eq. |7]is reasonably valid, as shown in Fig. |3] 

We point out that in the context of more general field theories of related time- 
evolving systems, apart from the AM, two other bosons observed in CDW systems 
are of interest. The first is created as a result of the annihilation of domain walls 
shown in Fig. 6. Upon annihilation, two finite-frequency (massive) ?/;— waves are 
emitted which propagate perpendicular to the surface into the sample and towards 
the surface. Calculations of the "wave effect" on the temporal evolution of the AM 
spectral shape appear to be confirmed by experiments in DyTea (discussed here) and 
inTbTes [12] • 

The second bosonic excitation of interest is the phonon mode (PM) with an 
equilibrium frequency of uipM =1.68 THz. The time-evolution of this mode is very 
different from the evolution of the AM: it does not show critical behaviour near t^- 
Instead its frequency is just slightly renormalized in the low symmetry state. Before 
tc, the PM frequency is ujpm =1-85 ± 0.02 THz. Because its frequency ujpm in the 
ordered state is lower than ujam (= 2.2 THz), it crosses the AM as the latter hardens 
after tc as indicated in Fig. 5. At Ati2 ~ 1.3 ps the AM interferes with the PM, and 
the phonon frequency is renormalized for Aii2 > 2 ps, to loam = 1-68 ± 0.02 THz. 

In experiments where the control parameter is temperature T rather than time, 
this behaviour is well understood jS]: below the critical temperature Tc, as the AM 
hardens with decreasing T, it crosses the PM at some intermediate temperature, 
displaying mode mixing. The coupling of the AM and the PM in compounds such as 
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DyTes is relatively weak, the off-diagonal matrix element being S ~ 0.1 ± 0.01 THz 
[S]. The PMs form a large reservoir of excitations which are weakly coupled to the 
order parameter. The symmetry of the PM plays an important role in the temporal 
behaviour. By definition, the order parameter is totally symmetric, (A representation) 
so even-symmetry modes are expected to couple strongly to it. When the system 
has inversion symmetry, such as here, (DyTea has D2h point group symmetry) the 
odd-symmetry modes do not couple to the tp to first order. Odd symmetry modes 
thus display some characteristics of weakly interacting massive particles in cosmology. 
The observed PM has a large mass (frequency) and is weakly coupled to the order 
parameter, either because they have inappropriate symmetry (in which case the weak 
coupling comes from higher order interactions), or the matrix element coupling to A 
symmetry excitations is small. In contrast to the AM and the fermionic excitations 
which emerge after the SBT takes place, the PM excitations exist before the SBT (i.e. 
before and after the Big Bang within the cosmological paradigm). By analogy with 
the PM, dark matter excitations may be thought of as a remnant excitations from 
before tc- As such they are external to the GL (or Standard) model. 

We conclude that building on the analogy between a laser-induced e-h plasma 
in condensed matter systems and the plasma created in the collision of elementary 
particles, or with primordial conditions in the Universe, gives us a practical laboratory- 
scale playground for the exploration of temporally evolving systems though symmetry 
breaking transitions. Remarkably, the CDW model system give rise to two kinds of 
cosmological analogies, namely of the AM as the analogue to the Higgs boson and odd- 
symmetry phonons as dark matter excitations. The CDW systems reveal some hitherto 
unexplored bosonic excitations created under nonequilibrium conditions which may be 
expected to have observable counterparts in other temporally evolving systems. The 
reflectivity response function derived for a superconductor allows the investigation 
of the trajectory of the order parameter through the transition, opening the way to 
studies of intrinsic vortex dynamics. 
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